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We present evidence of D"-D'^ mixing using a time-dependent amplitude analysis of the decay 
D" -> K+-K~Tr° in a data sample of 384 fb~^ collected with the BABAR detector at the PEP-II 
e+e^ collider at SLAC. Assuming CP conservation, we measure the mixing parameters s^^^o = 
[2.61 to.68(stat.) ± 0.39(syst.)l%, y^^^o = [-0.06 l^;g (stat.) ± 0.34 (syst.)]%. The confidence 



level for the data to be consistent with the no-mixing hypothesis is 0.1%, including systematic 
uncertainties. This result is inconsistent with the no-mixing hypothesis with a significance of 3.2 
standard deviations. We find no evidence of CP violation in mixing. 

PACS numbers: 13.25.Ft,11.30.Er,12.15.Mm,14.40.Lb 



The mixing between the flavor eigenstates \D'^) and 
\D^) of the neutral D meson depends upon the mass 
and width differences of the mass eigenstates. In the 
standard model (SM), the D°-D° mixing contribu- 
tion from short-distance effects is negligible pj. This is 
due to Glashow-Iliopoulos-Maiani suppression of the first 
two quark generations and Cabibbo-Kobayashi-Maskawa 
suppression of the third. Long-distance effects from in- 
termediate states that couple to both i)° and D° domi- 
nate. Their contributions to the mixing parameters are 
difficult to predict, but are estimated to be of the order 
10~'^-10~^ [l|. Several recent studies report evidence for 
mixing parameters at the 1% level [2]. This is consistent 
with some SM expectations and provides constraints on 
new physics models [3|. If mixing occurs, the physical 
eigenstates |£'i,2) = p\D^) ± 91^") must have different 
masses Mi^2 or widths ri.2. Here |pp -I- \q\'^ = 1. Mixing 
is parametrized using x = 2 (Mi — M2)/iTi + T2) and 
y = {Ti- T2)/{Ti + T2), where 1 (2) refers to the al- 
most CP-even (odd) eigenstate. If CP is conserved, then 
\p/q\ = 1 and arg(q/p • Af/Af) = 0. Here Af {Af) is the 
amplitude of the transition of the Z?° (-D°) to the final 
state /. 

In this letter, we analyze events in which the flavor 
of the Z)° [J] is measured at production. We present 
the first time-dependent amplitude analysis of the D^ — > 
K~^'k~it'^ Dalitz plot to extract the mixing parame- 
ters. Previously, we studied the time dependence of 
D^ -^ i^r+TT^vr" decays integrated over large regions of 
the Dalitz plot. We found no evidence for mixing [5|. 
However, certain regions of the phase space are more 
sensitive to mixing than others (see below). There- 
fore the amplitude analysis is potentially more sensitive 
than our previous work. The wrong-sign (WS) decays 
D^ -^ K^Ti^TT^ and the more copious right-sign (RS) 
decays 13° — ^ X^ 7r+ 7r° are reconstructed. The RS de- 
cays proceed through a Cabibbo-favored (CF) amplitude. 
The WS decays proceed through a coherent sum of a 
doubly Cabibbo-suppressed (DCS) amplitude and a CF 
amplitude produced by mixing. We identify RS and WS 
decays by reconstructing the !?*+ -^ D^ir'^ , D^ -^ Kittt'^ 
decay chain. The flavor of the _D° candidate is deter- 
mined from the charge of the low- momentum pion (tt^). 
The DCS and the CF amplitudes are described with iso- 
bar models [a] as outlined below. 

The time-dependent decay rate is a function of the 

Dalitz variables S12 = it^'^+tt- ^^^ '^13 ~ ^k+tt°' ^^ 
depends on the DCS amplitude Ajr{si2, S13) = {f\H\D^) 

and the CF amplitude y4j(si2,si3) = {f\H\D°) l], where 



/ = K'^TT tt". In the limit \x\, \y\ <^ 1, it is given by: 

dNjr{si2,Sl3,t) 



dsi2dsi^dt 



-^*{|^/f 



lAjll^djl [ycosSf -xiiin5f\ (Fi) + 
l±^\Ajnvtf] (1) 

where <5y^(si2, S13) = arg[^*-(si2, si3)Aj(si2, S13)]. The 
first term in Eq. [T] is the DCS contribution to the WS 
rate; the third term is a pure mixing contribution; the 
second term arises from the interference between DCS 
and mixing CF amplitudes. We determine the CF am- 
plitude ^jT in a time-independent Dalitz plot analysis 
of the RS decay sample, and use it in the analysis of 
the WS sample. The DCS amplitude Aj is extracted 
along with the mixing parameters. In the isobar ap- 
proach. At and At are described as a coherent sum of 
amplitudes, where each amplitude accounts for a reso- 
nance contribution. From inspection of the Dalitz plots 
(Fig. [T^-b), WS decays proceed primarily through the 
resonance D^ -^ A'*+7r^, while RS decays proceed pri- 
marily through D° -^ K^p+ [1]. For both Af and Aj, 
one complex amplitude must be fixed arbitrarily. 

The advantage of an amplitude analysis is that the 
interference term in Eq. [1] produces a variation in aver- 
age decay time as a function of position in the Dalitz 
plot that is sensitive to the complex amplitudes of the 
resonant isobars as well as the mixing parameters. In 
this study, the change in the average decay time and 
the interference between the Z3° — + i4r*+7r~ and D'^ -^ 
p~ K~^ amplitudes are the origin of our sensitivity to 
mixing. Our analysis is sensitive to x and y in the 
form y'l^^^o = y cos Sktv-^o - x sin Sktttt" and x'j^^^o = 
X cos (5/f ttttO + y sin Sxt^t^o where (5if Tru-" is the strong inter- 
action phase difference between the DCS D^ -^ p~K~^ 
and the CF Z?" — > K^p~^ amplitudes, and cannot be 
determined in this analysis (note that Sj^^^^o is in gen- 
eral different from the analogous phase in other decays). 
Mixing is implied through a non-zero value of x'^^^o or 
of y^^^o- We define ^/(si2, S13) = roA^'^^{si2, S13) and 

^/(si2,si3) =^j^(si2,si3)- Here 



rl = Nws/ {Nrs ■ / Nf{si2,si3,t)dsi2dsi3dt 



(2) 



where Nws (^fls) is the number of WS (RS) events in 
the sample. The variation of average decay time in the 
Dalitz plot due to mixing depends on the ratios x'^^^o /fQ 
and y'l^T^^o/fo'i these are the parameters that we extract 



directly from the data and use to determine the signifi- 
cance of the mixing resuh. 

The amphtudes entering the WS analysis are de- 
scribed as a sum of isobar components Aj that are 
parametrized with Breit-Wigner functions, A^r — 

Y.'jZi^'"^^ aje-""^' Aj{m]^+^_,m]^+^o), where aj and 5j 
are the strong interaction amplitudes and phases of the 
f^' resonant amplitude [6|. For the K-tt S-wave compo- 
nent we use a parametrization derived from K-ir scat- 
tering data [7|, which consists of a Xq(1430) resonance 
together with an effective non-resonant component. 

We analyze a data sample of 384 fb~^ collected with 
the BaBAR detector [1] at the PEP-II e+e~ colhder 
at SLAC near a center-of-mass energy of 10.58 GeV. 
Charged tracks are reconstructed with a silicon-strip de- 
tector (SVT) and a drift chamber (DCH), both in a 1.5 T 
magnetic field. Particle identification is based on mea- 
surements of ionization energy loss {AE /Ax) in the SVT 
and DCH together with measurements from a Cherenkov 
ring-imaging device. Photon energies are measured with 
a CsI(Tl) calorimeter. All selection criteria, the fit proce- 
dure and the systematic error analysis are finalized before 
we search for evidence of mixing in the data. 

Selection criteria are based partly on those of Ref. [51 
and are identical for the RS and WS samples. We re- 
quire the 7r+ candidates to have a transverse momentum 
pLAB ^ 0.12GeV/c, where LAB indicates the labora- 
tory frame. We reject electrons that mimic tt^ using 
dE/dx measurements. We use kinematic selection crite- 
ria to eliminate electrons from pair conversions. The en- 
ergies of photon candidates used to form tt*^ candidates 
are required to be greater than 0.1 GeV; the invariant 
mass of photon pairs forming a tt" must be in the range 
0.09 < TO^o < 0.16 GeV/c'^. We require the tt" momen- 
tum p^i^^ to be greater than 0.35 GeV/c. The recon- 
structed invariant mass for the D'^ candidates must have 
1.74 < TOk^^o < 1.98GeV/c2. The tt^ and D° masses 
are then set equal to their nominal values [9| and the 
D* is refitted 10] with the constraint that its production 
point lies within the beam spot region. The Z?*+ invari- 
ant mass and D^ measured decay time tp^^^o are derived 
from this fit. We require 0.139 < Am < 0.155 GeV/c^ 
where Am = mK-K-n'^-Ks ~ "^^K-mr"- To reject D* candi- 
dates from B decays, we require the D'^ center-of-mass 
momentum to be greater than 2.4 GeV/c. For events that 
contain multiple D* candidates with shared tracks, the 
candidate that yields the largest fit probability for the de- 
cay chain is retained. The three-dimensional flight path 
is used to determine tx-K-n" and its uncertainty ut- For 
signal events, the typical value oi at is 0.23 ps; we accept 
D* candidates with at < 0.50 ps. 

We extract the signal and background yields from a 
binned extended maximum likelihood fit to the mi^j^jro 
and Am distributions (Fig. [TJ:-d). For subsequent anal- 
ysis, we retain D* candidates in the signal region de- 



fined as 0.1449 < Am < 0.1459 GeV/c^ and 1.8495 < 
m^^^o < 1.8795 GeV/c^. Our final RS (WS) sample 
is composed of 658,986 (3009) events with a purity of 
99% (50%). The efficiency of the signal region selection 
is 54.6%. 





ml^, [GeV V] 



m|.^.[GeVV: 




m^.„„[GeV/c^] 



0.14 0.142 0.144 0.146 0.148 0.15 
,2n 



A m [GeV/ci 



FIG. 1: Dalitz plots for the (a) RS and (b) WS D'^ samples. 
The reconstructed (c) D" mass and (d) Am distributions for 
the WS sample requiring respectively (c) 0.1449 < Am < 
0.1459 GeV/c^ and (d) 1.8495 < m^^^o < 1.8795 GeV/c^ 
The fit results used to extract the yields are shown by the su- 
perimposed curves. The light histogram represents the mistag 
background, while the dark histogram shows the combinatoric 
background. 

The RS sample is used to determine the CF isobar 
model parameters af^ and df^ , as well as the decay 
time resolution function, which is parametrized as a sum 
of three Gaussian functions with a common mean, with 
widths given by the per-event at times a different scale 
factor for each Gaussian. The reconstructed RS D" sig- 
nal decay time distribution (Fig. [5^) is described by a 
probability density function (PDF) consisting of an ex- 
ponential function convolved with the resolution func- 
tion. The resolution function parameters and D'^ lifetime 
are determined in an unbinned maximum likelihood fit. 
The mean value of the resolution function is found to be 
4.2 ± 0.7 fs. This value is consistent with the magnitude 
expected from instrumental effects, and the associated 
systematic uncertainty is determined by setting the value 
to zero. As a cross-check we determine the D° mean life- 
time to be [409.9 ± 0.8 (stat. only)] fs, in agreement with 
the world average [410.1 ± 1.5 (stat. -f- syst.)] fs [9|. 

The D° candidates in the WS signal region can be di- 
vided into three categories: signal events, combinatorial 
background, and incorrectly tagged RS events (mistag), 
each one described by its own PDF, whose parameters 
are determined in an unbinned maximum likelihood fit. 



During the fit procedure, the number of events in each 
category is fixed to the value obtained from the fit to the 
mxn-TrO and Am distributions. 

The PDF describing the WS time-dependent Dalitz 
plot is given by Eq. [1] convolved with the ti^-^^o resolu- 
tion function. The at PDFs for signal and background 
are taken from the RS data. The DCS amplitudes and 
phases for each resonance and the mixing parameters 
are determined in the fit. The CF Dalitz plot ampli- 
tudes arising from mixing are taken from the fit to the 
RS sample described in the previous paragraph. The 
mistag events are parametrized using an empirical PDF 
obtained from the RS data, since mistag events contain 
correctly reconstructed RS D^ decays. The PDF describ- 
ing the combinatorial background is constructed by av- 
eraging the (si2, Si3, ij^^Tro) distributions obtained from 
the WS mxTTir" sidebands: this accounts for correlations 
between those three variables that might be present in 
the data. 

The results of the time-dependent fit of the WS data, 
the cij"~'^ , Sj"-^^ and fit fractions fj [6(, are given in 
Table HI The fit fraction of the non-resonant contribu- 
tion to the K-TT S-wave is absorbed into the ii'Q^(1430) 
and i^Q°(1430) fit fractions. Projections of the fit results 
are shown in Fig. [2)3-d. The change in log likelihood 
(—2 A hiC) between the fit with mixing and with no mix- 
ing i^K-Kno/'^o = y^^^o/ro = 0) is 13.5 units, including 
systematic uncertainties. For two degrees of freedom, 
the confidence level that the result is due to no-mixing is 
0.1%. The significance of the mixing result is equivalent 
to 3.2 standard deviations, and thus constitutes evidence 
for D'^-D'^ mixing. 



Using a Bayesian approach, by integrating the likelihood 
function with respect to x'j^^^o and y^u-Tr"' assuming 
a flat prior distribution, we obtain x'j^^^o =[2.6ltog8 
(stat.) ±0.39 (syst.)]% and y^^^o = [-0.06l°j^ (stat.) 
±0.34 (syst.)]% with a correlation of —0.75. 
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FIG. 2: (a) Proper time distribution for RS events with the fit 
result superimposed. The distribution of background events is 
shown by the shaded histogram, (b) Proper time distribution 
for WS events, (c, d) m^+^_. and ?ri^+^o projections with 
superimposed fit results (line). The light histogram represents 
the mistag background, while the dark histogram shows the 
combinatoric background; 



TABLE I: Fit results for the WS D° data sample. The total 
fit fraction is 102% and the x^ /ndof is 188/215. The resuHs 
for x'j^^^a /to and y'j^^^o /ro include statistical and systematic 
errors; their total linear correlation is —0.34. 



Resonance 



p(770) 

7^2*° (1430) 

i^*+(1430) 

7<'*+(892) 

ii'o°(1430) 

7^*° (892) 

p(1700) 



,UC'i' 



1 (fixed) 
0.088 ±0.017 

6.78 ± 1.00 
0.899 ± 0.005 

1.65 ±0.59 

0.398 ± 0.038 

5.4 ±1.6 



jr^ 



(degrees) 



(fixed) 

-17.2 ±12.9 

69.1 ± 10.9 

-171.0 ±5.9 

-44.4 ±18.5 

24.1 ±9.8 
157.4 ± 20.3 



/. (%) 



39.8 ±6.5 
2.0 ±0.7 
13.1 ±3.3 
35.6 ±5.5 
2.8 ±1.5 
6.5 ±1.4 
2.0 ±1.1 



x'j^^^a/ro = 0.353 ± 0.091 ± 0.052 
y'K^^a/ro = -0.002 ± 0.090 ± 0.059 



To derive the values of a^'^^^o and Vxtt-k" ^^ ^^^* deter- 



mine Tr 



^ [5.25 



+!^;3^ (stat.) ±0.12 (syst.)] 



1-3 



xlO "^ usmg 
Eq.m We then generate 10^ (a:^^^o/ro,y^^^o/ro) points 
in accordance with the fit covariance matrix, assuming 
Gaussian errors (width given by the total uncertainty 
including systematics) . For each point, we compute tq 



using Eq. [5] and determine values for 



'K-n 



and 



2/_ff7r7r"- 



Extensive validation of this fitting procedure is per- 
formed using Monte Carlo (MC) experiments based on 
the PDF shapes and DCS amplitudes extracted from 
data. The validation studies are performed over a wide 
range of mixing parameters. These studies demonstrate 
that the fit correctly determines the mixing parameters 
to within a small offset of 0.2-0.3cr, where a is the statis- 
tical uncertainty. These small biases are a consequence of 
the relatively small size of our data sample and become 
negligible if MC samples with higher statistics are used. 
We correct the final result for this offset. 

Sources of systematic uncertainty for 

^K-KTrol'^oiy'K-K-Kol'^o), related to the choice of the 
isobar model and the experimental assumptions, are 
considered. For each effect we refit the data with an 
alternative assumption and extract the overall correlated 
uncertainty for the fitted parameters. We estimate 
the Dalitz model uncertainties [0.38cr (0.35cr)], where 
a is the statistical uncertainty, by varying the mass 
and the width of each resonance within their error and 
by using alternative parametrizations for the isobar 
components Aj in the fit: the largest error arises from 
uncertainties in the K* and p parameters and from 



uncertainties in the parametrization of the K-tt S-wave. 
Systematic uncertainties related to the number of signal 
and background events [O.lScr (0.22a-)] are evaluated by 
varying them according to their statistical uncertainties. 
Similarly, the definition of the signal region, the at 
requirement, and the selection of the best D* candidate 
are varied. The effect on the mixing parameters is 0.50cr 
(0.37ct). Variations in efficiency across the Dalitz plot 
contribute systematic uncertainties of 0.09a (O.lOcr). 
The txTiTrO resolution function parameters are varied 
within their errors. The offset is also set to zero. The 
systematic effect is O.llcr {0.09a). The total systematic 



error on x 



K-ttttO 



ol'To iVKTr-Ko/'^o) is 0.57(7 {0.66a). 



The same procedure is applied separately to the WS 
ZJ^-tagged (+) and D°-tagged (— ) events to search for 
CP violation in mixing. We find x'^^^o = (2.53 
0.39)%, y'+^^o = ( 



-0.05;^-^? ± 0.50)%, X 



+0.54 I 
0.63 ^ 



-0.67 



KTrir" 



{3.55t'oil ± 0.65)% and y'^^^„ = (-0.54^?:^^ ± 0.41) 



-0.83 

respectively, and thus observe no evidence for CP viola- 
tion. The correlation between x'i „ {x'^ o ) and y'^ „ 
(2/m) is -0.69 (-0.66). 

In conclusion, our data are inconsistent with the no- 
mixing hypothesis with a significance of 3.2 standard 
deviations. Our results thus constitute evidence for 
mixing. For the mixing parameters we find x'j^^^q = 
(2.6ltl]:^^ ± 0.39)% and y^^^o = (-0.06t°j5 ± 0^4)% 
with a correlation of —0.75. These values are consistent 
with our previous result [5| and with some SM estimates 
for mixing. No evidence for CP violation is found. 
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